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Abstract 

The partition function and the one- and two-body distribution functions are evaluated 
for two hard spheres with different sizes constrained into a spherical pore. The equivalent 
problem for hard disks is addressed too. We establish a relation valid for any dimension 
between these partition functions, second virial coefficient for inhomogeneous systems in a 
spherical pore, and third virial coefficients for polydisperse hard spheres mixtures. Using 
the established relation we were able to evaluate the cluster integral 62 (V) related with the 
second virial coefficient for the Hard Disc system into a circular pore. Finally, we analyse the 
behaviour of the obtained expressions near the maximum density. 

1 Introduction 

The Hard Spheres (HS) and Hard Disks (HD) systems have attracted the interest of many physi- 
cists owing to they constitute prototypical simple fluids. Hard-core models, such as HD and HS 
are paradigmatic entropy driven systems. The search of their accurate equation of state is a 
long-standing problem of great importance in statistical mechanics [HElGilll]. The evaluation 
of virial coefficients [TJ \S[ O \7\ [H [9] and the analysis of phase transitions [H [10] are some of the 
principal subjects. The monodisperse systems as well as bidisperse and polydisperse ones have 
been extensively studied. 

Recently special attention was paid to systems of few HS and HD confined in small vessels. The 
study of these inhomogeneous systems has shed light on aspects of loss of ergodicity [Tl] , freezing 
and glass transitions [12| . thermodynamic second law [13], and other fundamental questions of 
statistical mechanics and thermodynamic [14l 115]. 

Exact and approximate analytic properties of hard spherical bodies at different dimensions 
have been very important in the evolution of Free Energy Density Functional Theories for HS 
[TBI [17] . Besides its own relevance the HS system is the starting point of several approximate 
theories of liquids [Hill], therefore exact analytical results become still more interesting. Even 
though the apparent simplicity of HS and HD fluids, a few exact analytical results are known at 
present. For homogeneous systems the third and fourth virial coefficients for HS and HD have been 
calculated on the monodisperse EDI ED E21 E3] and polydisperse systems EHESIES], numerical 
results for higher order terms have also been obtained EH EH]- In addition, for inhomogeneous 
systems of HS into a spherical cavity the analytic expression for the second virial coefficient is 
known [29] . Until present the exact canonical partition function of few bodies in a pore has been 
solved for two HD [30J and recently for two HS [31], in a rectangular box. Such results for equal 
sized particles were compared with dynamical simulation approaches [15, 32, 33). The three HD 
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in a rectangular box problem was analyzed [34| but the complete integration of its configuration 
integral (CI) could not be performed. Exact analytic results exist for the zero dimensional limit, 
which corresponds to a so narrow cavity that is able to contain one particle at most. Also, several 
exact results have been obtained for the Hard Rod system (HR) the ID version of HS [U[35] which 
actually may be considered completely solved. 

The aim of the present work is to evaluate the canonical partition function for two spherical 
particles in a spherical pore at two and three dimensions. Expressions for the free energy and 
one particle density are also derived. In Sec. [2] of this work we show that in the canonical 
ensemble a spherical pore containing HS or HD can be seen as another particle. We establish 
a relation between the CI of the canonical partition function for N polydisperse hard spherical 
particles and the respective CI for N-l particles inside a spherical vessel. In Sec. [3]we evaluate the 
canonical partition and density distribution functions for two HD and HS into a spherical pore. 
The obtained expressions apply to the non additive system of unequal sized particles in a pore, 
which includes the more usual and restrictive additive system. The analysis of the CI expressions 
and the thermodynamic properties of the system are shown in Sec. 21 Finally in Sec. [5] we present 
our conclusions. 



2 The In- Out Diagram relation 

In this section we will refer to non additive HS, but the discussion also applies to HD and to the 
equivalent system in any dimension. We are interested in a system with fixed number of particles 
(N) and temperature, then we will consider the canonical ensemble. In such ensemble the partition 
function factorizes into CI and a trivial kinetic terms. The CI of the system is characterized by the 
set of hard repulsion distances between each pair of particles {d^}. The system is additive only 
if the repulsion distances are dij = Ri + Rj for positive radii of particles {Ri}. When the above 
relation is not fullfilled or any other particular assumption is made on the set {dij} the system is 
non additive. We will analyze the CI for a 7V-HS system in an open (and virtually infinite) volume 
(Qn)- In order to obtain a non null contribution to the integral, each pair of particles ij must 
fulfil Tij > d^, where r.y is the distance between particles ij. The above relation between and 
dij is a consequence of the exponential Boltzmann's factor. Then we have 

Qn = I ■■■ J Yl e ij dr i--- dr N , (1) 

<ij> 

e- — pxp( — d) -(r -)) — 1° if ^ < ^ (overlap) 

e„ -ex P ( <p v [r v )) ^ if > ( no over | ap ) 5 W 

where 4>ij(rij) is the hard core spherical potential for the pair of particles ij. In virial series and 
related subjects this factor is usually named /, however in the framework of canonical partition 
function it appears naturally as an exponential function. Then we prefer to name it e or 
function or simply an e-bond, which will be plotted as a dashed segment. Such & Qn can be 
expressed as a Ree-Hoover's graph (also called modified star graph) [HI [20] 
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where the N-th particle at Eq. J3]) was drawn separately and its e-bond appears explicitly. The 
remaining (N-l )-particle system with e-bonds between pairs is outlined as a rectangular box. The 
relation between the Mayer / function (or /-bond) and e is given by 

eij = 1 + fij , (4) 



Qn = \ i 

> i 
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fij — 



— 1 if rij < dij (overlap) 
if rij > dij (no overlap) 



(5) 



We can transform each e-bond that links the pairs Ni (i : 1, . 
relation Eq. |(4|) performing the decomposition over particle N 



, N — 1) through the usual bond 
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In this representation all the equivalent graphs on the left of Eq. that takes into account per- 
mutations of the N-l particles were omitted for clarity. It is possible to map this general problem 
to the 7V-HS additive system defined by the set of particle radii {Ri} through the assignment 
Ri + Rj. We define the additive/subtractive generalization of the additive system by the 



\Ri ± Rj\, and specifying which set of pairs assume the minus sign. If we are 



assignment d, 

interested in a system of additive spheres contained in a spherical vessel, as is our case, the ex- 
clusion distances between the center of the pore (which is the N-particle) and the center of each 
i-particle is duq = (P — Ri), where P is the radio of the spherical pore. We name to this specific 
additive/subtractive {dij} as an additive-in-pore system. 

Through the assignment of {dij} for an additive-in-pore system, the last graph of Eq. ^ may 
be interpreted as the CI of (N-l) additive HS in a pore (Qm V -_ 1 0, due to the condition nN < djjv 
with 1 < i < N 



1 imposed by the /-bonds in Eq. ©. Therefore {day} in that graph should 

p 

(N-l)' 



be named inclusion (and not exclusion) distances. More precisely, last graph is times Q p 



being Voo the volume of the infinite space, 

Q(N-i) = / ■•• / U<iN>fiN Yl <i:j> e l jdri...dr N - l dr N , 

= fp'-'J I\ <lJ> e lJ dr 1 ...dr N - 1 , 
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where the integration domains depend on pore size P. From now on, last graph of Eq. J6]) will 
be called the (TV-i)-in-pore graph, even in the more general framework of non additive systems. 
If we have N — 1 equal (and additive) hard spheres of radio a/2 into a pore of radio P, on the 
additive-in-pore picture we assign duq = (P — a/2). The same system can be analyzed as an 
entirely additive system with c?ijv = [P — o) + a/2, if P > a. We can obtain a Mayer type 
decomposition by transforming iteratively each e-bond in Eqs. (J3j [Tj) by an /-bond. 



3 Two bodies in a pore 

In this section we will apply the in-out relation for three bodies. In Eq. we perform the simple 
decomposition over the particle P and show the relation between the 2-particles-in-pore graph 
and the Mayer type graph 

® © 

©— ® ©— ® 



(8) 



All the integrals in Eq. ^ may be directly evaluated, the easiest are the simply connected 
graphs that become factorable [36J . In the evaluation of different graphs we will use the notation 
of [26]. The exclusion/inclusion distances between pair of spheres are AB, PA, PB. Being 
that {AB, PA, PB} are completely independent we deal with a non additive system. Now we 
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Z = Z(r, Ri,R 2 ) 


Z{(r, R U R 2 ) 


Z' n {r, Ri,R 2 ) 


Z"(r, R U R 2 ) 


HS 




Sp(Ri) - Z 


S P (R 2 ) - Z 


Voo - SpiRi) - S P (R 2 ) + Z 


HD 




Cr(fii) - Z 


Cr{R 2 ) - Z 


V 00 -Cr{R 1 )-Cr{R 2 ) + Z 



Table 1: Relation between the set of zeta functions {Z, Z^,Z[j,Z"} and the integration of one 
dark particle linked with two other white particles. Note that in Z[ and Z' n , R\ is the repulsion 
distance with the left white sphere, while R 2 corresponds to the right one, with R\ > R 2 . 

will briefly analyze the two bodies or one body in a pore problem. Assuming that PA is the 
exclusion/inclusion distance between A and P, we have for HS 

© ® = Sp(PA) , (9) 

&—®=V 00 -Sp(PA), (10) 

A-jr 

Sp(R) = -jR 3 , (11) 

where Sp(R) is the volume of the sphere of radio R. The Eq. |(9|) is the accessible volume for a 
particle in a pore with inclusion distance PA (eventually PA = P — A), which is related to the 
second coefficient of the pressure virial series [26] for a mixture of two HS gases. The CI of two 
spheres with a repulsion distance PA (eventually PA = P + A) is given by Eq. lflO|) , where Voo is 
the volume of the total space. These obviously are not novel results. 

3.1 Two spheres in a pore 

Now we focus on the two spheres in a pore problem and then we need to integrate the 2-in-pore 
graph (left-hand side graph on second row of Eq. ((HJ). Before integral evaluation we may state 
that some set of values {AB, PB, PA} produce a zero sized graph, since A and B spheres may 
not enter into the pore P. If AB > PA + PB, then Eqs. (f5]) and |(5]) imply that tab > tpa +tpb 
thus avoiding the triangular relation and therefore the 2-in-pore graph becomes zero. From now 
on we will assume AB < PB + PA and PB > PA, which is compatible with a pore of radius P 
that contains two particles of radius A and B (P > A + B) with A > B. 

For the calculus of density distribution it is useful introducing the function Z{r, R\, R 2 ) |26J, 
which represents the volume of intersection of two spheres of radius R\ and R 2 separated by a 
distance r (assuming R\> R 2 ) 

( Sp{R 2 ) iovr<R 1 -R 2 
Z(r,R lt R 2 ) = \ I 2 (r,R 1 ,R 2 ) other (12) 
[ for r > i?i + R 2 , 

h(r, Rx,R 2 ) - ^-{Ri +R 2 - rf{r 2 - 3{R 1 - R 2 f + 2r {R r + R 2 )) . (13) 
12 r 

The function I 2 (r, Ri,R 2 ) is the volume in the partially overlapping configuration, and is symmet- 
ric through the permutation Ri <-> R 2 . Other functions related with Z(r, R\,R 2 ) are summarized 
in Table [TJ (the full extended form of Table Q] appears in Appendix [A} . This set of functions is 
related with the integration of one spherical body linked with other two. As we have two types 
of bonds e and / there are four different Zeta functions. The use of these functions accounts us 
to write down two relevant distribution functions in a simple way: the pair distribution function 
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g2(r AB ) (in which the position of the pore center was integrated) and the one body distribution 
function of one particle pi(r A ). Both, r A and r B have their origin in the center of the pore 

92{r AB ) = (QZy 1 e AB Z(r AB ,PB,PA) , (14) 
P \-i, f \ / Z{(r A ,PB,AB) for AB < PB 



= m- i-f PA ) x | ^ a ,ab7b) ^ ab ; pb , ( 15 ) 

f ^ (nP\~ 1 ( t wJ ^(ra.PAAB) for AB < PA 
Pi(r B ) = (Q a ) (-/pb) x ) ' (16) 



Z' u (r B ,AB, PA) for > PA , 

where pi(r B ) is equivalent to pi(r J 4) through the transformation A «-> B in Eq. l[T5j) . Performing 
the complete integration (for example, in Eq. (flA)) ). we found the CI of two HS in a pore 

pPA+PB 

Qa= / 47r r 2 ^ Z (r AB , -PB, PA) dr AB , (17) 

JAB 
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(PA + PB - AB) 6 (AB 3 + 3AB 2 (PA + PB) 



P J -(PA + PP + 3AB)((PA-PP) 2 -2PAPB)) for AB > PB - PA 
i%^PA 3 ( PB 3 - AB 3 ) for AB < PB - PA. 



i2 



(18) 
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Thus, using Eq. lf20j) we may evaluate the CI of two additive spheres in a pore 
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which is valid as long as AB < PA + PB, whereas = if AB > PA + PB. There are explicit 
relations between the non additive, the additive, and the additive-in-pore systems. These are 
expressed in Eqs. f(T9|) and i[2"0j) respectively, which allow the map from one to other system when 
the triangular relation is fulfilled 

1 + B 2A = AB + PA-PB 

(19) 



(20) 



Q2 = (yf) {P - A- Bf (P 3 + ZPAB - A 3 - B 3 ) , (21) 

above expression is valid for P > A + B, whereas Q P = if P < A + B. With the aim of check, 
using Eqs. ([H [TBI [H]) we have obtained the third virial coefficient in the additive poly disperse 
system [25} [26] the right-hand side graph on second row of Eq. (JSj) . In addition, in the case of 
equal sized spheres with radii A we obtain 

QP=f^\ (P-2A) 3 (P 3 + '3A 2 P-2A 3 ) , (22) 

Q P 2 = (f ) 2 (P - Af {^j (2A) 3 (P - Af + ^ 2 (2A) 4 (P - A) 2 - ^(2A) 6 . (23) 

These two equations are equivalent. In Eq. ([22[) it seems clear the existence of a root on Q P 
at P = 2A. On the other hand, Eq. ([23| shows the volume and area dependence, being V = 
(4tt/3)(P — A) 3 and Ar — An(P — A) 2 . A simple relation connects Q P with the second cluster 
integral 62(F) (which was stated in the second row of Eq. ((Sj) ) ; due to 62(F) has been analytically 
evaluated in previous works [29] we were able to check the validity of our result. Distribution 
functions for two additive spheres in a pore are shown in Appendix [Bj 
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3.2 Two disks in a pore 



The two dimensional problem of two HD in a circular pore is very similar to the problem of two 
HS in a spherical pore in three dimensions. Therefore we will briefly outline the aforementioned 
analysis concerning HD. In the two particle system or one disk in a pore, we obtain 



© © = Cr(PA) , 



0— ® = V« 
Cr(R) = 



- Cr(PA) 
ttR 2 , 



(24) 

(25) 
(26) 



where Cr(R) is the surface area of the circle with radio R. To maintain an unified point of view 
for both systems, HD and HS, we name the measure of the total space Voc (although it is actually 
an area) . The Eq. (f24|) is the accessible area for a particle in a pore with an inclusion distance PA 
(eventually PA = P — A), which is related to the second coefficient of the pressure virial series 
on a mixture [26]. The configuration integral of two disks with repulsion distance PA (eventually 
PA = P + A) is given by Eq. ([25)1 . For the two disk in a pore problem we have 



Z(r, R U R 2 ) 



Cr(R 2 ) 
h(r, Ri,R 2 ) 




for r <Ri-R 2 
other 

for r > Ri + R 2 



I 2 {r,Rx,R 2 ) = R\ arccos f + 2 f Ri Rl ] + R\ arccos f ~ P }^ 1 \ 

~\ y/{r +RT+ R 2 ){-r + R x + R 2 ){r - R[ + R 2 ){r + R[ - R 2 ) 



(27) 



(28) 



The function I 2 (r,R\,R 2 ) is the surface area in the partially overlapping configuration. The set 
of functions related with Z(r, R\, R 2 ) for disks are summarized in Table Q] (see the full extended 
form in Appendix [Aj. On terms of the set of Z functions it is possible to write down the pair 
distribution function g 2 (r AB ), and the one body distribution function pi(r A ) 



9i(tab ) = {Q 2 ) e AB Z{r ABl PB.PA) 



Pi(r A ) = {Q P 2 y 1 {-fp A ) x j |j 



Z[(r A ,PB, AB) for AB < PB 
(r A ,AB, PB) for AB > PB 



(29) 
(30) 



where Pi{tb) is equal to pi(r A ) under the transformation A «-» B in Eq. f30|) . Both, r A and r B , 
have their origin in the center of the pore. Performing the complete integration of Eq. (J29j) , for 
example, we find the partition function for two disks in a pore 



Q p 2 



PA+PB 



2irr AB Z(r AB ,PB, PA) dr AB , 



(31) 



AB 



Q p = { 



ttk (PB 2 + PA 2 + AB 2 ) 
+7T PB 2 PA 2 (tt - arcco S [ PB2 +ff A AB2 }) 
— 7f PB 2 AB 2 arccos[- PB2 - pA2+AB ^ 
l arccos\- 



2 PB AB J 
-PB 2 +PA I +AB 2 
2PAAB 



tt 2 PA 2 (2PB 2 - AB 2 ) 



for AB > PB 
for AB < PB 



PA 
PA 



(32) 



k = -y/(PB + PA + AB) (-PB + PA + AB) (PB — PA + AB) (PB + PA — AB) 



(33) 



where K represents the area of the triangle with sides PB, PA, AB. The arguments of the arccosine 
function in Eq. lj32j) are the internal angles of this triangle and their sum is n. In particular, for 
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the additive in pore system with different sized disks, the previous equation can be simplified 



Q P = 2ttk(P 2 +A 2 +B 2 -P{A + B) + AB) 

+AP Bf{P - A) 2 (n - arccos[ p2 -;}^_-* B ]) 
-*(P - B) 2 (A + Brarccos[ p f-_ B ^l + 4 B ] 
-*(P - Af{A + B)*arcco8[ - p $Z$i£2 ) AB ] , 



(34) 



k = \j P B A(P - A - B) . (35) 

The Eq. {34J) applies in the case P > A + B, whereas Q P = if P < A + B. Using Eqs. <f8j fl9j 
and [34l) the third virial coefficient is correctly obtained in the additive polydisperse system [24J. 
By evaluating the above expressions for two disks with equal radii, we find 



arccosl — — 



or 



Q P = 2tt A^JP {P- 2A) ((P - A) 2 + 2A 2 ) 

+ 2tt(P-A) 2 ((P-A) 2 -(2A) 2 j, ; ,,, j j 

Q P = (±§V2irA 3 / 2 )(P-2Af 2 

x (l + ^j(P-2A) + ^ I (P-2A) 2 + 3 (P-2A) 



(36) 



(37) 



Eq. |(36]) has a root at P — 2A. Thus, in Eq. (J37J) we study the local behaviour of Eq. (J36J) near 
its root. Distribution functions for two additive disks in a pore is discused in Appendix [51 



4 Results 

Let us make a few remarks on the new results obtained in this work. The configuration integral 
in Eqs. (fl8ll32"|) are analytical functions of the system parameters in the domain AB, PA, PB > 
except when we focus on the non additive regime for AB = PB — PA. From here on we will 
restrict ourselves to AB < PA + PB. For HS the configuration integral is polynomial, although 
it is not true for HD. In the two additive bodies in a pore system, the CI Eqs. 11211 l34j) become 
null at close packing pore size Pep = A + B, even though systems do not lose the rigid rotation 
degrees of freedom. Over all density range P € {Pcp , +oo) CIs are monotonic analytic functions, 
therefore, none ergodic/non-ergodic transition happens and no van der Waals loop is expected 

[Sol M- 

As we are interested in thermodynamic properties we analyze the Free energy of the system. 
Taken into account the kinetic factor, which contains the temperature dependence and usual 
thermodynamic relations we obtain the entropy S 

S/k = -f3F + f3U = ln(Q P ) , (38) 

where U is the energy of the ideal gas (j3U = 2, 3 for two HD and HS, respectively). An important 
result derived from Eqs. l|2"H l3"4 |) is that the Free Energy diverges logarithmically at Pcp- In order 
to obtain the pressure we need to define the system volume V sys , both magnitudes are related by 



f3P w 



dln(Q F ) 



dln(Q F ) 



dP 



(39) 





where V^ ys is the area of the surface enclosing V sys . A particular volume definition implicitly 
induce a surface definition and then only modifies the surface area in Eq. 1(39]). Although we are 
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considering a small and very simple system composed by solely two hard spherical particles in a 
pore, it is not obvious which is the system volume. In principle, three different volumes may be 
considered, the volume accessible to the particle A, that is Va; the same for the B particle i.e. Vb 
(see Eqs. (|9] 124]) ); and the volume of the free space in the pore. The later may be interpreted as 
the accessible volume for a particle with a vanishing radio and does not depend on the size of the 
chosen particle. We propose that Va and Vb are the relevant volumes. In such a case the pressure 
on the wall became 



p fw = 



+ dln(Q£) 

0,V B dVB 



(40) 

0,Va 



dV A 

0P w = Pl {rA = P-A)+ Pl (rB=P-B) , (41) 

where Eqs. l|40] HTjl are generalizations of Dalton's Law, and the equivalence between them is 
trivial. 

For the sake of clarity from here on we will analyze the system of two equal sized particles, i.e. 
B = A and Pop = 2A An interesting point about Eq. (|23"]) for two HS, is the fact that it has the 
simplest possible shape. Indeed, each term has a very simple mean. The first term refers to the 
ideal gas, while the second term is clearly related to the first virial series correction (due to the 
interaction between particles) both of them appear in the homogeneous two body system. The 
third term is the product between the area of the system and the second surface virial coefficient 
at a planar hard wall W2 |37[ [38] , Finally the fourth term forces the root of CI at P = Pcp- 
From Eqs. ((221) and ([37]) the CI goes to zero at P C p with multiplicity 3 and 5/2 for HS and HD, 
respectively, and its value for two HD in a square box is 4 |30j. The mean of such values are 
probably related to the number of degrees of freedom lost at the ultimate solid or cp density. The 
logarithmic divergence of the Free Energy and Eq. (|40]l induce an order one pole at the maximum 
density. The equations of state for HS and HD are, respectively, 

and 

^=27^((^ + l^( P - 2 ) + °^- 2 ■ ™ 

Here, we choose A as the unit length. Note that Eqs. i|42]l and |43]) are truncated forms of the 
equation of state but the exact complete expressions can be easily derived from Eqs. 1(22] [36 ] [40]) . 
It seems that some of the afore analyzed properties about the ultimate solid are valid for any 
finite number of particles N. That is, the CI must go to zero at some P — Pcp(N) value with a 
multiplicity at-i(N), which is a positive real number. This characteristic implies the logarithmic 
divergence of the Free Energy and induces an order one pole in Pw at the maximum density. Then 
we conjecture that the equation of state near the ultimate solid density, takes the form 

^ = n D (p- ip-i ( (F^b) + "° + ai{P ~ Pcp) + ° 2{P - Pcp) ) ' (44) 

where D is the dimensionality of the system and is the solid angle integral. This last equation 
can be rewritten in the following way 

PPw^ = wn(i-(Vo/vy/ D y 1 + a ^ D -^ 

(45) 

+ ao +MPcP-i) (F _ Pcp) + 02{p _ Pcp) 5 

where V — tin (P — ^) D /D and V a = fin {Pcp ~ I) 13 '/D. It is interesting to mention that first 
term of the Eq. lj45]) is equal to the equation of state for the bulk solid phase of HS and HD 
systems at the highest density limit |39J. 

As was mentioned before, CIs are monotonic analytic functions, then the pressure is a mono- 
tonic analytic function of P € (2, +00) and does not develop a van der Waals loop. This is not the 
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P/A 



Figure 1: (Color online) Equation of state for two HD and HS. Continuous curves correspond 
to HD (up) and HS (down) calculated in the present work. For comparison, we also include the 
results for: two HD in a square box |30j (dashed line), and Monte Carlo simulation for two HS in 
a spherical pore [40] (circles). 

case for two and three HD in a rectangular box |15 [ [30 ] [34] . In Fig. [T]we display the equations 
of state 11421 143|) as a function of pore size, both of them diverge at P = 2. For comparison, the 
equation of state for two HD system confined in a square box pore is also plotted, there the half 
side of the square box pore was taken as the size parameter P. The equation of state for two 
HD in a square box has two branches that are joined on a non analytic point at P = 2, while 
the divergence appears at P — Pep = 1 + 1/V2 and takes place when the discs are fitted on the 
two opposite vertices of the box. These two equations of state for two HD on different shaped 
pores converge as the pore size grows. Present results of f3P\v for HS perfectly agree with the 
Monte Carlo results [40] (extracted from Fig. 1 on that work) which were obtained by numerical 
integration of the same equations. 

The density distribution of one HS particle corresponding to a system with two equal sized HS 
(see Eq. f56|) in the Appendix [Bj has been plotted on Fig. O Density profiles for several radii 
of the pore P are shown in this figure by a solid line, each curve has a maximum at r = P — 1 
and fall off to zero discontinuously for r > P — 1 . Density distributions change gradually from the 
complete radial localization, when P — Pep = 2 to the quasi uniform distribution at P — 4. For 
a pore size P > 3 the density at r = becomes non null and a constant density plateau develops 
for < r < P — 3. The dashed and dot-dashed curves show the evolution of the minimum (at 
r = 0) and the maximum (at r = P — 1) values of the density profiles as a function of pore radio 
P. The dotted curve indicates once again the contact density p(P,r = P — 1), which is related 
to the pressure of the system by Eq. (|4T|) . but now as a function of P — 1. These three lines are 
functions of pore size, and refer to the top abscissa axis. We find that all these properties are 
also exhibited by the two HD in a circular pore system. We can mention that at P = 3.4069 the 
central density reaches its maximum value, which is close but greater than half of the density at 
contact. 

The complete dependence of the second cluster integral 62 (V) for a system of particles contained 
into a spherical vessel is related to the second virial coefficient by B(V) = —b2(V). From Eq. 
we obtain 



where V = (P - A) 3 , b 2 = -^-(2A) 3 , and Ar — Air (P — A) 2 . As was discussed before (in 




(46) 



Using Eq. (|22| we derived 62 (V) for two equal HS in a spherical pore, 




(47) 
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P/A i 




r/A 

Figure 2: (Color online) Density profiles for two HS system as a function of r for various P 
values. Both, P and r are in units of A. Continuous lines indicate results for pore radii P = 
2.01, 2.5, 3, 3.5, 4. Other curves are referred to top abscissa axis. Dashed curve shows density 
at central point p(P, r = 0) , whereas dot-dashed curve corresponds to density at contact p(P, r = 
P — 1). The dotted curve is the shifted dot-dashed one (see text). 

Sec. I3.ip this expression was previously calculated in literature [29]. For the first time, we present 
the complete dependence of b 2 (V) for a HD system into a circular vessel. Then, using Eqs. 1|36|1 
and |46| 

Vb 2 (V) = b 2 V+l(2AfAr~±n(2Af^^ + 3 (- F ^j , (48) 

where V = tt (P - A) 2 , b 2 = -\k{2A) 2 , and Ar = 2n(P-A). Actually V is the accessible 
area (see Eq. (|26l) ) and Ar is its perimeter. Although in Eq. ({48}) we choose express b 2 (V) as a 
truncated series expansion of the complete expression, the exact form of b 2 (V) (for two HD and 
HS with equal or unequal size) may be obtained with the same procedure. The factor multiplying 
Ar is the second virial coefficient of the surface tension (in fact the linear tension) W 2 [37} [38], as 
far as we know it was not previously evaluated. The third term is the first finite size-curvature 
correction and it has a minus one power on [P — A), which is two units smaller than the former 
term as happens in HS at Eq. l|4*T)) . All the forthcoming terms of the series (not shown in Eq. 
(HE}) are negatives with odd degree in (P — A)~ . The expansion of Eq. lj48j) converges so quickly 
that, if we truncate at fifth order, the deviation is smaller than one percent in the worst case for 
P = 2A, when the system attains its highest density. 

We can mention that CI of two Hard Rod system confined into a segment (the one dimensional 
equivalent of two HD and HS in a spherical pore) is easy to evaluate analytically [21 [35]. Then, 
this work provides a set of analytic properties for dimensions D = 1, 2 and 3, and any general 
approximate theory for inhomogeneous confined liquid phase [17] may be compared with these 
exact results. 

5 Conclusions 

We have established a deep relation between CI of a N-polydisperse system of hard spheres, CI of 
a (N-l)-polydisperse system in a spherical shaped pore, and Mayer type diagrams |19[ 136]. The 
relation valid in any dimension was used in the present work only for checking purposes. A similar 
relation was stated previously in the Gran Canonical Ensemble [29j. 

The canonical partition and distribution functions of two HS and two HD in a spherical shaped 
pore were analyzed in an exact analytic framework. The analysis presented here was made assum- 
ing a very general system of two HD and HS, with different sizes and non additive potential. The 
obtained exact CI expression allow us to evaluate the thermodynamic observables of the system 
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(Free Energy, Energy, Entropy, Pressure, etc) [H1I36]. In addition, we investigate the analytic 
properties of CI and the equation of state at the ultimate solid or cp density, and we compare our 
results with the known exact result of two HD confined in a square box. Based on this study we 
present a general proposal for the equation of state at the ultimate solid, in a system of N hard 
spherical particles contained into a pore. We also find the second virial coefficient for a HD system 
in a circular pore. Future works will involve the analytical evaluation of the CI for two particles in 
other simple shaped pores, and the analytic study of the CI for three particles in spherical pore. 
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A Appendix A: Expanded definition of Z[, Z[ Y and Z" 

In this appendix it is showed the expanded form of the Z(r, i?i, R2) related set of functions for 
spherical bodies. We name D to the dimensionality of space and we assume that R\ > R2. In 
the present notation an unified point of view for spheres and disks is chosen. The family of Z 
functions in arbitrary dimensions reads, 

for r<R!-R 2 

Z(r.Ih-lh) { („/?,../?,) other (49) 

for r > Ri + R 2 , 

for r<R 1 ~R 2 

Z,'(r../,V/? 2 ) = { J\'"iJh)-iy ) (r,R 1 ,R 2 ) other (50) 

for r > Ri + R 2 , 

for r<R 1 -R 2 

Z' [l [r.n l ..ll 2 ) = { l\ l "ill.)-I^ D) (r,Ri,R 2 ) other (51) 

for r > i?x + R 2 , 

Ko - (Ri) for r<R 1 -R 2 

Z"(r,fli,Jfc) = { V 00 -l[ D) (R 1 )-l[ D) (R 2 )+I { 2 D \r,R 1 ,R 2 ) other (52) 
Voo - l{ D) {R 2 ) - l[ D) (R 2 ) for r > R 1 + R 2 , 

where i[ D \R\) is the volume of the D-sphere with radio R\ (the expression for D — 3 and D = 2 
are in Eqs. (|HI26p ). and I 2 D \r, R\,R 2 ) is the volume of intersection of two D-spheres with radii 
Ri and R 2 , separated by a distance r in the partial overlapping configuration (the expression for 
D = 3 and D = 2 are in Eqs. (|13|28p ). 




B Appendix B: Distribution functions in the additive-in-pore 

system 

Assuming an additive-in-pore system, i.e. AB = A + B, PA = P — A and PB = P — B, the 
distribution functions for two HD or HS into a spherical pore (from Eqs. (fl4l [15]) and |29l l30|) ) 
take the form 

ao(r^)- (On' 1 x I I { 2 D) (r AB ,P-B,P-A) if A + B < r AB < IP - (A + B) 
92 (r AB )-{Q 2 ) X I Q i£r AB <A + Boir AB >2P-(A + B), 

(53) 
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max[l[ D) {P- B) - l[ D) (A + B),0] if < r A < abs[P - A-2B] 
Pi{r A ) = (Q2) lx { l[ D) {P - B)- I ( 2 D) (r A ,P- B,A + B) if abs[P- A-2B] < r A <P-A 
" if r A > P - A , 

(54) 

the P\{r B ) function can be found by the permutation of symbols A <-> B in Eq lj54j) . If we have 
equal radii spherical particles, expressions become 

ao(r*A-(O p V l xl 4 D \r AA ,P-A,P-A) if 2A < r AA < 2(P — A) , , 

92{r AA ) - { Q 2 ) x| q . f ^ < ^ or r ^ > 2(p _ A) ^ (55) 

{maa;[7 1 (D) (P - A) - l[ D) (2A), 0] if < r A < abs[P - 3A] 
l[ D) {P- A)- 4 D) (r A ,P- A,2A) if abs[P-3A] < r A < P - A 
" if r A > P - A . 

(56) 

Although the existence of the exact one body distribution function for two HS in a spherical pore 
is mentioned in [41], neither explicit expression was shown, nor further analysis was performed 
there. 
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